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ABSTRACT 

The method devised by Lewis et al. (1993) for calculating the light curve of a mi- 
crolensed point source is expanded to two dimensions to enable the calculation of light 
curves of extended sources. This method is significantly faster than the ray shooting 
method that has been used in the past. The increased efficiency is used to obtain much 
higher resolution light curves over increased time scales. We investigate the signatures 
arising from different source geometries in a realistic microlensing model. We show 
that a large fraction of high magnification events (HMEs) in image A of Q2237+0305 
involve only one caustic, and could therefore yield information on the structure of the 
quasar continuum through the recognition of a characteristic event shape. In addition, 
the cataloguing of HMEs into morphological type will, in theory, enable the direction 
of the transverse motion, as well as the source size to be obtained from long term 
monitoring. 

Key words: gravitational lensing - microlensing - quasars - numerical methods. 



1 INTRODUCTION 

Much of the study of cosmological microlensing has focussed 
on the modelling of the images of Q2237+0305 in which mi- 
crolensing was first observed by Irwin et al. (1989) and Corri- 
gan et al. (1991). Many authors (e.g. Paczynski 1986; Wamb- 
sganss & Paczynski 1991; Witt, Kayser & Refsdal 1993; 
Lewis et al. 1996) have used different numerical methods 
to model variations in the observed flux of this object. Sta- 
tistical arguments have then been used to place constraints 
on the size of the emission regions of the quasar. The sim- 
ulations have also shown how long term monitoring of mi- 
crolensing events can yield information on the composition 
of the lensing galaxy. As an example, Schmidt & Wambs- 
ganss (1998) place limits on the lower end of the mass range 
from an observed lack of microlensing in Q0957+561. 

Interpretation of model light curves for extended 
sources calculated using ray shooting (e.g. Kayser, Refsdal & 
Stabell 1986; Wambsganss & Paczynski 1991) has been hin- 
dered by the poor statistical coverage of the simulated light 
curves. The majority of magnification maps have been com- 
posed of 500x500 or 1000x1000 pixels. Therefore when the 
resolution is high enough, the light curves have only covered 
a few months or years for source sizes of < 10 16 cm, instead of 
the centuries required to adequately sample high magnifica- 
tion events. In addition to the ray-shooting method, approx- 
imations to microlensed light curves for extended sources 



have been obtained in two ways: through combining the 
parametric representation of caustics with the ray-shooting 
method (Wambsganss, Witt & Schneider 1992), as well as 
through integrating an extended profile over a line of point 
source magnifications (Witt & Mao 1994) . The first of these 
methods computes the light curve at a caustic by convolv- 
ing the near caustic approximation (Eqn ^) with the source 
intensity profile. It therefore assumes the validity of Eqn ^ 
for a complex caustic network. The second approach does 
not consider the behaviour of the magnification pattern in 
the direction perpendicular to the source line, and will thus 
overestimate the proportion of HMEs which involve only one 
caustic. 

High magnification events can be classified into a num- 
ber of different classes. If the event is produced by the source 
crossing a single caustic (SHME), then the characteristic 
shape of the caustic will contain information on the geome- 
try of the quasar emission regions (Fluke & Webster 1998). 
On the other hand, high magnification events produced by 
the source crossing two or more fold caustics (MHMEs) are 
harder to interpret. Classification of HMEs, and the subse- 
quent physical interpretation, requires full modelling of in- 
dividual events. This places limits on the coarseness of the 
microlensing simulations that may be of use. The observed 
microlensed event in image A of Q2237+0305 in 1990 had 
a rise time ~ 26 days. Therefore as a rough guide, a tem- 
poral resolution of ~ 1 day is required for the successful 
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modelling of a well sampled HME in Q2237+0305. Lewis 
et al. (1993) and Witt (1993) independently introduced a 
method for calculating the images of a line, and therefore 
the light curve of a point source, microlensed by a collec- 
tion of point masses. Section ^| describes how this method 
can be expanded to calculate the microlensed light curve 
of an extended source with an arbitrary surface brightness 
profile. In section ^| this extended method is used to com- 
pute microlensed light curves of simple model quasars with 
different intensity profiles. The number and type of high 
magnification events are then determined. Throughout the 
paper, standard notation for gravitational microlensing is 
used. The Einstein radius of a 1M© star in the lens plane is 
denoted by £ ; in the source plane it is r\ \ the normalised 
shear due to mass in the galaxy that is outside the region 
important for microlensing is 7 and the convergence or op- 
tical depth is k. Where required, a cosmology having Q — 1 
and A = with H = 75 km sec -1 Mpc -1 is assumed. 



2 CONTOURING ALGORITHM FOR 
EXTENDED SOURCES 

2.1 Construction Of The Light Curve For An 
Extended Source 

The method used to compute the light curve of an extended 
source is based on the contouring method. A series of mi- 
crolensed light curves are computed. These light curves have 
parallel source tracks that together cover the source and pro- 
duce a regular grid of points in the source plane. A source 
intensity profile is then be convolved with this grid to pro- 
duce a total observed magnification of the source. When 
an extended source is microlensed, its magnification at the 
point 77 is calculated using the convolution: 



J p(ff' -jKgHdY 

J P(fj' -fj)d?rf 



(1) 



where P(ff' — if) is the intensity profile of the source centred 
on fj at the position f}' , and fJbpirf) is the 2-dimensional mag- 
nification pattern. If the magnification pattern is calculated 
using the ray shooting technique, then the computed mag- 
nification at the coordinates of a pixel is the average mag- 
nification over the region of source plane covered by that 
pixel. Thus the convolution ([!]) can be easily evaluated nu- 
merically as a weighted average. However, if the contouring 
method is used to build a magnification grid in the source 
plane, then the magnification at a given grid point is that 
of a point source. The singularities in f-i p (rf) at 



0}, 



(2) 



where £ is the image position, must be integrated out di- 
rectly as the magnification of the point source is not nec- 
essarily a representative average of the region surrounding 
that point. This is the most important aspect of the calcu- 
lation as much of the flux and therefore the information in 
the light curve originates from the parts of the source which 
lie very close to a caustic. 

To calculate the two dimensional integral (Eqn |l|) over 
the finite grid, the 1-dimensional integrals for each line of 
the grid are first evaluated. The two dimensional integral 



is then evaluated by integrating these areas perpendicularly 
to the direction of the source line. The number of lines that 
cover the source strip was determined as follows. The trans- 
verse velocity is unknown, however as the event rise time was 
measured at ~ 26 days, we know that ~ 50 points per source 
width will sample to the same level as potential observation. 
We therefore set 50 lines as the minimum coverage for the 
strip. 

The convolution of the source intensity profile with each 
source line is performed in sections. A source which lies on 
the side of a caustic line such that it has two associated 
critical images is referred to as being inside the caustic. The 
parts of the line that do not lie just inside a caustic are 
smooth and so can be integrated in a straight forward man- 
ner. For sections immediately inside a caustic the magnifi- 
cation is highly nonlinear with distance. To overcome this 
nonlinearity we compute the numerical integration across 
the caustics with the help of the analytical result of Chang 
& Refsdal (1979) and Chang (1984). The magnification p, v 
of a point source at a small perpendicular distance d from a 
fold caustic is 



(3) 



where \i is the contribution from non-critical images, 0(d) 
the Heaviside step function, and a is a constant . This result 
assumes that locally the caustic is a straight line. It also as- 
sumes that the source position remains approximately con- 
stant with respect to any other nearby caustics so that the 
value of fj, is also constant. In the calculations that follow, 
Eqn|^ is used at distances < 10~ 2 t; o , that can be compared 
to the curvature of the caustics and the inter-caustic dis- 
tances of ~ lr] a . Fluke & Webster (1998) have obtained a 
higher order approximation to Eqn fa that treats the shape 
of the caustic as a parabola. A possible direction for future 
work will be to use this result together with information on 
the caustic curvature to obtain a better approximation to 
the magnification during both the interpolation and inte- 
gration algorithms. 

To utilise Eqn ^ in the evaluation of Eqn |l] we need to 
know the location of the intersection of the caustic and the 
source lines. This is determined during the contouring algo- 
rithm, when the approximation Eqn ^| with fi = is used for 
interpolation of magnifications between two source points in 
the high magnification regime (Lewis et al. 1993). We im- 
plement the contouring algorithm as described by Lewis et 
al (1993). The image solutions are found by stepping along 
the image curves in the one direction. Therefore the case of 
the image line crossing a critical curve is determined from a 
sign change in the magnification corresponding to a parity 
flip of the image. This sign change results in the direction 
of travel of the corresponding source line solution being re- 
versed at the caustic, the location of which is then found 
through a linear interpolation in /ip 2 . During the contouring 
algorithm the intersections of all caustics with each source 
line are recorded for use in the convolution (Eqn ^]) . 

Figure ^ shows a schematic representation of the source 
grid in relation to the caustic. The labels of source positions 
and magnifications which follow refer to this figure and de- 
scribe their relationship to the source grid. The determined 
positions of the caustics y ca ust are used in conjunction with 
the values of y"a,y~b and flp(ya),IJ^,(yb) a ^ the two nearest grid 
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Figure 1. Left: Representation of of the points on a finite source grid, with the labelling as referred to in the text, and Right: 
Representation of a curve fitted to the amplifications at grid points on the middle source line which crosses the caustic structure shown. 
The positions of the grid points are labelled by the light dashed lines. 



points (inside the caustic), and the straight caustic approx- 
imation (Eqn ^) to obtain the values of fi , and a for the 
current source point, expressions for which are 

„ _ fJ-piVa) - HpiVb) ,., 

— -, — 1 — _ 1 \ > 

(\Vcaust — Va\)~2 — (\y ca ust ~ 3 

and 

Ho = Hp(Va) ~ — — T - (5) 

[\ycaust -y a \) 2 

Note that the value of a obtained is not identical to that in 
Eqn ^| but is weighted by a factor of \/l/ sin9 where 9 is the 
angle between the source trajectory and the tangent to the 
caustic. This factor cancels with that relating the perpen- 
dicular distance d with the distance to the caustic along the 
source line (\y C aust — y a \)- Since the singularities in fJ-piv) 
at the caustics are removed in the initial source line inte- 
grations, this approach is clearly only valid for integrations 
along source trajectories that cross caustics having no part 
within the source strip tangential to the source line. Before 
integrating out the singularity, the convolution routine must 
locate the two grid points straddling the caustic, and decide 
which of the two points is inside the caustic and which is 
outside. 

The presence of a caustic between two adjacent points 
(h,l) and (h±l,l) on a given grid line is determined using 
the following rules; firstly, the magnification /x p (h,l) of the 
current point must be a local maximum that is greater than 
Hp — 10, and secondly, the ratio between the point source 
magnifications at the current grid point and that immedi- 
ately in the forward /u p (h+l,l) or backward /u p (h— 1,1) di- 
rection must be greater than 10. This second condition en- 
sures that the integration of Eqn [I] using the approxima- 
tion (Eqn |^) is not carried out in a high magnification but 
smoothly varying region such as that located outside of a 
cusp. In addition, the second condition determines whether 
the current source point is immediately inside or immedi- 
ately outside the caustic, and so which grid points should 
be used in the evaluation of Eqns H and pi The situa- 



tion of non-isolated caustics is recognised by the algorithm 
when the difference between the position of the current grid 
point and the recorded position of the corresponding caustic 
(\yi a — yi caust\) is greater than the separation between adja- 
cent grid points (yidiff = \yia — yib\) on the source line. In 
this case the interpolation scheme (Eqn ^) has failed during 
the contouring algorithm. The routine also looks for another 
caustic lying along the source line within 3 grid points of the 
first, and if either of these are the case, discounts the con- 
tribution to the flux from the point closest to the current 
caustic since the integration of Eqn |l] will not be valid. It 
should be noted that the routine was observed to be robust 
against variations in the parameters used in these tests. 

2.2 A Comparison Between Light Curves 

Computed Using the Ray-Shooting and 
Contouring Techniques. 

The backward ray tracing or ray shooting technique of 
Kayser, Refsdal & Stabell (1986) provides an excellent test 
of the validity of a light curve computed using the contour- 
ing method. An example of a microlensed light curve for 
an extended source, calculated using contouring technique 
(solid line) is shown in Fig |j| Also shown is the correspond- 
ing curve calculated from ray shooting (dotted line), and 
the numerical difference between the two (lower panel). The 
corresponding caustic network is shown in the top panel to- 
gether with the extremities of the source strip and allows 
features in the light curves to be related to locations on the 
caustic network. Details of the simulation parameters are 
described below. 

The light curves are calculated for the approximation to 
the bolometric luminosity of a thin accretion disc with outer 
radius R s — 0.0577 o (see Grieger, Kayser & Refsdal 1987). 
The surface brightness profile is described by the following 
expression, where r is the radius: 

Paccr(r) = f^-) 3 for ^ < r < R s . 
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Figure 2. Middle: The light curve from the accretion disc source calculated using the contouring (solid line) and ray shooting (dotted 
line) algorithms. The convolution (Eqn |l|) was performed on a 50 line grid, the source scale is in units of r\ and the horizontal dashed 
line is the average magnification for a model with parameters re = 0.2, 7 = —0.2. Bottom: The difference between the two light curves. 
Top: The corresponding caustic network. The extremities of the source strip are plotted as dashed lines superimposed on the caustic 
structure. 







otherwise 



(6) 



This profile was chosen as a possible representation of the 
profile of the continuum region in an active galactic nucleus. 
The source strips were computed with 50 source lines sepa- 
rated by 0.0027/ o , so that the convolution (Eqn|l|) for the con- 
toured curve was computed on a 50x50 grid at each source 
point ff. For the calculation of the ray-shooting source strip, 
the rays were shot in a regular grid in the lens plane, and 
the number of rays used was determined so that there were 
100 unlensed rays per source strip pixel. The region of lens 
plane where images were looked for was the same in both 
the contouring and ray-shooting calculations. The lens used 
in this case was 30 point masses with an optical depth of 
k ~0.2, and a normalised shear of 7 =-0.2 applied parallel 



to the source line. Calculations of these light curves using 
source grids composed of 100 and 200 lines were also made. 
Note that this simulation is for comparison between the ray 
shooting and contouring methods and is not intended as a 
realistic simulation at the quoted microlensing parameters. 
However the theoretical average magnification is shown su- 
perimposed on the light curves in figure ^ (dashed line) for 
comparison. 

The lower panel of figure |^ shows excellent agreement 
between the two methods of calculation, particularly at 
points where the source is not covering the intersection of 
two caustics or a cusp. In either of these situations, the ray 
shooting calculation is expected to be more reliable as it is 
not dependent on the use of the straight caustic approxi- 
mation. This approximation is not applicable at either the 
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Table 1. This table refers to the light curves shown in figure Q 
The ranges shown describe positions along the source track, and 
are positioned over the HMEs. The table shows the root-mean- 
square of the difference between the light curves for the accretion 
disc source when calculated with the contouring technique and 
with the ray shooting technique. 



Range of 
Source Line (^) 


rms 
(50 lines) 


rms 
(100 lines) 


rms 
(200 lines) 


-1.1 -0.9 


0.011 


0.004 


0.002 


-0.8 -0.6 


0.009 


0.005 


0.003 


-0.6 -0.2 


0.014 


0.008 


0.005 


0.0 0.2 


0.014 


0.007 


0.004 


0.4 0.6 


0.029 


0.018 


0.007 


0.65 0.8 


0.013 


0.007 


0.004 



intersection of two caustics or at source points just inside 
a cusp where two caustics are within three grid points of 
each other. In both instances the integration routine will fail 
since the interpolation scheme is based on this approxima- 
tion, and so the point source magnifications obtained from 
the contouring algorithm will be incorrect. As a quantita- 
tive measure of the deviation between the two methods of 
calculation, the root-mean-squares of the difference between 
the two light curves during high magnification events on the 
source trajectory shown in figure ^ were computed for the 
50, 100 and 200 line grids. These are listed in table [l] and 
show root-mean-square differences which are of order 1% of 
the size of the high magnification event for the calculations 
using a 50 line grid. The difference between the two methods 
of calculation is converging to zero for finer grid sizes. 

HMEs in microlensed light curves which involve only 
one caustic exhibit variations in their characteristic shape 
due to different source geometries (Fluke & Webster 1998; 
Schneider, Ehlers & Falco 1993) . This difference between the 
light curves for the two sources is not a function of the grid 
size, and we note that it is significantly larger than the nu- 
merical difference between the ray shooting and contouring 
methods. In addition we note that while the numerical un- 
certainty (ray-shooting minus contouring) is predominantly 
present in regions of the light curve where the gradient is 
largest, the difference in the shape of the light curves is 
present for the duration of the HME. 

The contouring algorithm has two major advantages 
over the ray shooting method. Firstly, it solves the lens 
mapping function, rather than finding the magnifications 
through a statistical calculation. Secondly, and more impor- 
tantly, there is an increase in speed of computation. For 
calculations of light curves from different source geometries, 
the magnification grid in the source plane must have a small 
resolution so that the convolution (Eqn |l|) can be computed 
over many grid points. To obtain an accurate magnifica- 
tion in a pixel using the ray shooting algorithm requires 
at least 100 un-lensed rays per pixel (Schneider, Ehlers & 
Falco 1993). Computation of a fine grid over a long source 
line therefore becomes impractical for large star-fields and 
small source sizes. Computation times on a Unix worksta- 
tion (Sun UltraSPARC 1) for the 50, 100, and 200 line source 
grids used to calculate the light curves in figure ^ were found 
to be ~ 50 seconds, ~ 100 seconds and ~ 200 seconds for the 



contouring algorithm, compared with ~ 150 minutes, ~ 600 
minutes and ~ 2400 minutes respectively for the same cal- 
culation using ray shooting with the deflection of rays being 
calculated across the minimum region required to collect the 
total flux (for this lens configuration). During realistic mod- 
elling the addition of a large number of point masses means 
that this region must be enlarged. This is discussed further 
in section 3.1. The difference in required CPU time is due 



to the number of times that each algorithm must compute 
the bending angle. The length of time required for the cal- 
culation of a given light curve is approximately independent 
of source size for contouring but is proportional to the in- 
verse of the square of the source diameter for ray shooting. 
In addition, computation time increases linearly with the 
number of source lines for contouring, but with the square of 
the number of source lines for ray shooting. The ray shooting 
method can easily produce a number of light curves from the 
same model starfield. However, the source tracks will have 
significantly overlapping shooting regions and will therefore 
not be statistically independent. Thus, relatively speaking, 
contouring becomes more efficient for smaller sources and 
finer grids. 

Haugan (1996) suggested the use of contouring to de- 
fine the regions of the source plane in which rays need to 
be traced so that they are only collected by the region of 
interest in the source plane, thus reducing the number of 
rays required. He obtains an expression for the efficiency of 
the hybrid method with respect to traditional ray shoot- 
ing in terms of the number of rays used (not including the 
contouring calculations). For a simulation of optical depth 
k — 0.2 this expression yields an efficiency of ~ 200 for the 
source track we have computed for this section. This number 
is comparable to the efficiency we obtained from computa- 
tion times for the 50 line strip. We note however that as the 
hybrid method still employs the ray shooting technique, its 
efficiency with respect to contouring is still proportional to 
the inverse of the number of lines in the strip, as well as to 
the source size. Thus contouring will still be more efficient 
for smaller sources and finer grids. 



3 MODEL LIGHT CURVES FOR Q2237+0305 

The substantial increase in computational efficiency allows 
a range of microlensing simulations to be explored. This en- 
ables the determination of statistical measures for compari- 
son with the data. 



3.1 Models with Static Stellar Positions 

The observation of microlensing in Q2237+0305 has 
spawned many attempts to characterise the variation in 
flux using numerical microlensing models. The ray shooting 
method has been used to produce caustic maps and light 
curves from various source sizes of each of the images in 
Q2237+0305. Examples of these are given in Wambsganss, 
Paczynski & Katz (1989). We base our models upon those 
presented in that paper. The models consist of point masses 
with positions distributed randomly in a disc, and masses 
determined from a Salpeter mass function p(m) oc m -2,35 
in the range Q.IMq < rrii < I.OMq. The model does not 
include a continuous matter distribution. The region of the 
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Figure 3. Sample light curves for the approximation to an accretion disc source (dark line) and a circular source with limb darkening 
(light line), together with the difference between the two curves. The microlensing parameters are k = 0.36 and 7 = +0.4. 



lens plane in which image solutions need to be found to en- 
sure that 99% of the total macro-image flux is recovered 
from all points on the source line is known as the shooting 
region. This region is found by first calculating the image 
of the source line, as lensed by a smooth screen of matter 
with the same optical depth as the starfield. The number 
of stars in the region about any point which collects 99% 
of the macro-image flux was calculated by Katz, Balbus & 
Paczynski (1986) and is given by 



N, = 100 



<m 2 > 



(m> 2 |(1-k) 2 - 7 2 |' 



(7) 



In the case of an applied shear these stars are distributed 
in an ellipse with a major to minor axis ratio of (1 — ft + 
7)/(l — K — 7). The shooting region is then given by the 
union of these ellipses centred on all parts of the image line. 
The minimum number of stars required in the model are 
contained in a disc having a radius R t which covers this 
shooting region. In the cases of 7 = —0.4 and 7 = +0.4 
we have chosen the radius of the disc of point masses to be 
1.2 x R, and i?* respectively. 

Initially, the stars are assumed to be stationary. The 
effect of including stellar proper mot ions in the microlens- 
ing model is discussed in section 3.3 The macrolens model 
of Schmidt, Webster & Lewis (1998) was used to determine 
the microlensing model, giving the parameters of k = 0.36 
and I7I = 0.40 as a description of the microlensing environ- 
ment of image A. A transverse velocity in the lens plane 
of 600 km s -1 was assumed (eg. Wambsganss, Paczynski & 
Schneider 1990) . The size of the simulations was then deter- 
mined such that the length of each was ~ 70 years for the 
model with positive shear, and ~ 140 years for the model 
with negative shear. Therefore the length of each of our sim- 
ulations is Wrjo where 7 = +0.4 and 20rj o where 7 = —0.4. 



The numbers of stars in our models are then 2032 in the 
case where 7 = +0.4 and 880 where 7 = —0.4. 

For Q2237+0305, the Einstein radius of a 1M© star in 
the lensing galaxy (2 = 0.0394) projected to the quasar red- 
shift (z — 1.695) is rj « 1 x 10 17 cm. Light curves were calcu- 
lated for the cases of 7 = —0.4 and 7 = +0.4 that correspond 
to source trajectories directed towards, and at right angles 
to, the centre of the lensing galaxy respectively. Contouring 
simulations were calculated in each of these cases for source 
radii of 5 x 10 15 cm(i? s ~ §a) an d 1 x 10 15 cm(.R a ~ jfa) 
using as the source profile an accretion disc model of the 
source (Eqn ^) , and an approximation to a circular source 
with limb darkening (Grieger, Kayser & Refsdal 1987), 



Pcir (?*) 







or < r < R s 
otherwise 



(8) 



The two source geometries have different surface areas but 
are normalised to have the same total unlensed flux. Again, 
the source grids were composed of 50 lines each. The light 
curves therefore contain 5000, and 25000 points for the 
^ a ~ 20 anc ^ ^ s ~ urn res P ec tively in simulations with 
7 = +0.4, and twice this number in simulations where 
7 = —0.4. This was considered a sufficient resolution to 
determine the presence of single caustic high magnification 
events, and other miscellaneous events due to clusters of 
fold caustics. For comparison of the two sources in a realis- 
tic model, examples of the light curves for both the R s — 2a 
accretion disc (dark line) and circular sources (light line) 
for the case of 7 = +0.4 are shown in figure [| The differ- 
ence between the light curves for the circular and accretion 
disc profiles is also plotted, and shows that differences are 
present in the model up to a level of 0.1 magnitudes. 
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Magnitude Change 



Figure 4. Point source magnification distributions for the pa- 
rameters k = 0.36 and 7 = +0.4 (dark line), 7 = —0.4 (light 
line). 

A random starfield will be homogeneous when consid- 
ered on a large enough scale so that the granular nature 
of its boundary will have a negligible effect. Wambsganss, 
Paczynski & Katz (1989), and others have used tree-codes 
and associated methods to calculate the bending angle due 
to fields of as many as 10 5 point masses. These codes have 
the effect of smoothing the matter far from the image so 
that only the contributions to the bending angle from the 
few hundred point masses nearest to the image position are 
summed directly. The magnification patterns produced in 
microlensing simulations respond in a complicated, and of- 
ten radical manner to small changes in the position or mass 
of one of the model stars (Wambsganss 1992). However, this 
change only occurs in parts of the source plane where the 
corresponding shooting region includes the point mass in 
question. A model starfield is therefore sufficiently large if 
it covers the shooting regions corresponding to all points on 
the source line. In the case of the models used in the current 
work, the simulations were sufficiently small that we found 
the use of tree codes to be unnecessary. 

To check that the parameters used in our simulations 
are appropriate we generated a sample of 100 point source 
light curves, in each of the cases 7 ±0.4 using the source line 
lengths and parameters described previously . We find that 
that the mean magnifications for the 7 > and 7 < light 
curves are 3.92 and 3.90 respectively. These values are con- 
sistent with the theoretical value of fx av = 1/((1 — k) 2 — r y 2 ) = 
4.00. In addition, the difference between theoretical and nu- 
merical values is consistent with that quoted for comparable 
models in (Lewis & Irwin 1995). This demonstrates that our 
models are indeed finding the micro images that account for 
the macro image flux. Figure ^ shows the probability that 
a point source in our models is magnified by a magnitude 
between m and m + 5m, where m — —2.5 x logio(fj,/ fj, av ) 
and fi is the microlensing magnification generated from our 
models in the cases of 7 > (dark line) and 7 < (light 
line). As required these curves are identical. In addition we 
can compare the distributions to those in Lewis & Irwin 
(1995). Our distributions are quantitatively consistent with 
those presented in that paper. 



3.2 High Magnification Events 

In order to analyse the high magnification events (HMEs), 
three different classes are defined. 

(i) Single-Caustic High Magnification Events (SMHE): 
HMEs that have the same general shape as the correspond- 
ing convolution of the profile with the straight caustic ap- 
proximation (see for example Schneider, Ehlers & Falco 
(1993)). The event must be complete, having a discontinuity 
in gradient at its leading edge as well as a drop in magnitude 
at its trailing end. It must include the two discontinuities in 
gradient, the small secondary peak, and the 4 changes in the 
sign of the second derivative of the curve. 

(ii) Multi-Caustic High Magnification Events (MHME): 
HMEs that have both a rise and a fall, as well as at least one 
discontinuity in their gradient. These are not recognisable as 
SHMEs as the source straddles two or more fold caustics. 

(iii) Cusp Events: HMEs where the source moves through 
a high magnification region outside a cusp. They have both 
a rise and a fall but no discontinuities in their gradient and 
are conspicuous in the model light curves because the shape 
of the HME is approximately independent of the source ge- 
ometry, and approximately symmetric about its maximum. 

The HMEs were classified interactively. 

Table ^ shows the average rates for the different types 
of HME in each of the four cases described earlier. These 
rates depend linearly on the transverse velocity. The errors 
in these, and results of later sections, were calculated as 
the Poisson error of the total number of events in the light 
curves. As was noted by Witt, Kayser & Refsdal (1993), and 
others, the rate of variability is greater for a source which has 
a trajectory parallel to the shear vector (7 > 0) due to the 
stretching of the caustic network. This change in rate is due 
to the fact that a source moving perpendicularly to the shear 
vector (7 < 0) will on average spend more time in regions 
of the source plane which contain no caustic clusters. For 
image A of Q2237+0305, and a small source the simulations 
of the current work show this change in rate to be a factor 
of approximately 2, in agreement with the results of Lewis 
& Irwin (1996). 

Since the direction of the source trajectory has a sig- 
nificant effect on the frequency of observed HMEs, there is 
a degeneracy between the effective transverse speed of the 
source and the direction of the source trajectory. In the case 
of Q2237+0305, it may be possible to eventually deduce the 
direction of motion relative to the caustic network by con- 
tinued monitoring of the four images, which are orthogo- 
nally positioned with respect to the centre of the lensing 
galaxy (eg. Witt & Mao 1994). However, the simulations 
show that we also expect that the ratio of the frequency 
of Cusp HMEs to total HMEs to be approximately twice 
as large for a source moving along a path perpendicular to 
the shear vector than parallel to it. Since this ratio is inde- 
pendent of the transverse speed, its measurement provides 
a direct method of determining the direction of motion. In 
addition this ratio appears to be independent of source size. 
The higher proportion of cusp events in the case of negative 
shear is expected because the preferential direction of for- 
mation of the caustics with respect to shear vector directly 
implies that the cusps also have a preferred direction. Iso- 
lated cusps form mostly on the outside of the long caustic 
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Table 2. The frequency of different classes of HME in the various classes of simulation. All measurements of time assume a lens plane 
transverse velocity of 600 km sec . 



Class of Total length of Cusp HMEs SHMEs MHMEs total HMEs fraction 

Simulations simulations (years) per decade per decade per decade per decade of SHMEs 



Rs = 


■ 05tj o 7 = 


+0.4 


1610 


0.12 ± .02 


1.21 ± 


07 


0.87+ .06 


2.20 + 


10 


0.55 ± 


.04 


R s = 


.Oljjo 7 = 


+0.4 


420 


0.18 ± .05 


2.47 + 


20 


0.38 + .08 


3.03 ± 


22 


0.82 ± 


.09 


R s = 


.05r? o 7 = 


-0.4 


1680 


0.13 ± .02 


0.63 ± 


05 


0.62 + .05 


1.38 ± 


08 


0.46 ± 


.05 


R s = 


.Olr/o 7 = 


-0.4 


700 


0.16 ± .04 


1.16 ± 


10 


0.19 + .04 


1.51 ± 


12 


0.77 ± 


.09 



clusters so that a source trajectory that is parallel to the di- 
rection of the caustic clustering will encounter more of these 
smooth high magnification regions. 

The total frequency of HMEs is a function of source size, 
with a higher frequency for smaller sources. This increase 
is not large for the range of source sizes considered here, 
and is due to the fact that a smaller source is less likely 
than a larger source to cover multiple caustics. In contrast, 
the relative proportion of SHMEs and MHMEs is a very 
sensitive function of source size, with SHMEs being ~ 6 
times as prevalent (with respect to MHMEs) for the R 3 = 

source than for the R s = ^ source. This result appears 
to be approximately independent of trajectory direction. 

In order to determine the statistical frequency of HMEs, 
frequent monitoring for a long period of time is required. 
The most important scientific result is the large percentage 
of HMEs in the model light curves that are also SHMEs. If 
the central region of the source quasar has a geometrically 
simple intensity profile then we would expect to see a char- 
acteristic shape in the observed light curves. This character- 
istic event shape could then be deconvolved to reconstruct 
the source intensity profile either by trial and error, or by 
a reconstruction method such as that suggested by Grieger, 
Kayser & Refsdal (1987). Information on the source intensity 
profile might be obtained after observing just a few HMEs. 



3.3 The Effect Of Proper Motion Of Stars 

Most of the numerical work on microlensing in Q2237+0305 
has concentrated on models of static starfields, with the flux 
variation resulting from an assumed effective transverse ve- 
locity in the lens plane of 600 km s -1 . The main motivation 
for this assumption is undoubtedly the enormous increase in 
computation time required to compute a light curve with the 
inclusion of stellar proper motions. This increase is the result 
of the entire magnification pattern having to be re-computed 
for each point on the model light curve. In addition to the 
positions of the point masses, the macro parameters are in 
principle also time dependent. However, these time varia- 
tions are negligible when compared with the motion of stars 
because changes in k and 7 occur over a scale of more than 
1000 years (Kundic, Witt & Chang 1993). 

The observed one dimensional velocity dispersion of the 
bulge in the galaxy of Q2237+0305 is « 215 kms -1 (Foltz 
et al. 1992). If this value is representative of the disper- 
sion at the position of image A, then it may have a sig- 
nificant impact on the statistics of HMEs. Two approaches 
have been taken towards the computation of the effect of a 
velocity dispersion on microlensing statistics. The velocity 



of the caustic network resulting from stellar velocities has 
been calculated by Schramm et al. (1993), and by Kundic, 
Witt & Chang (1993). The latter consider the area swept 
out per unit time by the caustic network produced from an 
ensemble of point masses with a velocity dispersion. They 
find that a one dimensional velocity dispersion will produce 
more HMEs than a static model with the equivalent trans- 
verse velocity. A second approach has been taken by Kundic 
& Wambsganss (1993) and Wambsganss & Kundic (1995). 
They investigate the microlensing contribution by comput- 
ing the change in the magnification pattern that occurs when 
the model stars are ascribed motions according to a velocity 
dispersion. Kundic & Wambsganss (1993) find a contribu- 
tion to HMEs from the velocity dispersion which is con- 
sistent with that of Kundic et al. (1993). For image A in 
Q2273+0305 Kundic & Wambsganss (1993) find that a one 
dimensional velocity dispersion of a ~ 215 km sec' 1 will in- 
crease the number of observed HMEs by about 10% if the 
transverse velocity is assumed to be vt = 600 km sec' 1 . 

Since stellar proper motions may have a large effect on 
the statistics of HMEs, it is necessary to investigate their 
effect on the existence and relative proportions of the differ- 
ent types of HMEs. Our prime motivation is to estimate the 
expected prevalence of characteristic SHMEs in an observed 
light curve. Therefore, rather than compute many points 
along a light curve, each with a sligh tly evolved starfield, 



many of the simulations from section 3.1 were re-run with 
the stars having been displaced in both the j/i and 1/2 di- 
rections by an amount determined from a speed distributed 
according to a Gaussian having a dispersion of 215 kms -1 , 
(note this assumes an isotropic velocity distribution), and a 
time equal to that required for a caustic with the projected 
bulk velocity to move l^x the source width. If a charac- 
teristic SHME is present in both the light curves, then it 
will also be present in a fully computed curve using a model 
with an evolving field. The time scale of this event would 
be either increased or decreased depending on whether the 
fold caustic was moving with or against the bulk velocity. 
The proportion of total HMEs that are SHMEs obtained in 
this manner will approximate the ratio SHME:HME for the 
model used. 

Table |^ gives the relative frequency of the different 
classes of HMEs when stellar proper motions are included in 
the microlensing model. These rates are obtained using the 
effective transverse velocity and do not take into account the 
increased frequency of HMEs due to the inclusion of stellar 
proper motions. The incorporation of proper motion has had 
the effect of reducing the proportion of SHMEs in the model 
light curves. This is a consequence of the evolving caustic 
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Table 3. The frequency of different HMEs in microlensing models that include stellar proper motions. All measurements of time assume 
a lens plane transverse velocity of 600 km sec" 1 . 



Class of Total length of Cusp HMEs in SHMEs in Total HMEs fraction 

Simulations simulations (years) both curves both curves per decade of SHMEs 

per decade per decade 



Rs — 


.Obrio 7 = 


+0.4 


770 


0.09 + 


.03 


0.89 ± 


.09 


2.26 + .15 


0.39 ± .05 


Rs = 


.Ol77o 7 = 


+0.4 


350 


0.24 + 


.06 


2.18 ± 


.20 


3.06 + .24 


0.71 ± .09 


R s = 


.05ri o 7 = 


-0.4 


1260 


0.09 + 


.03 


0.49 ± 


.06 


1.43 + .09 


0.34 ± .05 


R s = 


.01?7o 7 = 


-0.4 


700 


0.19 + 


.05 


0.94 + 


.11 


1.44+ .14 


0.65 ± .10 



network, however the frequency of SHME s is still an appre- 
ciable fraction of that obtained in section 3.1. Therefore the 
effect of proper motion of stars in the model is to increase 
the overall number of HMEs, but to decrease the fraction of 
SHMEs. 

Whilst the effective transverse velocity is unknown, and 
may well be relatively small, the fact that stellar proper 
motions can produce microlensing events at an appreciable 
rate, means that even in the absence of a transverse velocity, 
both HMEs and SHMEs are common. The presence of a 
SHME in one of the time separated light curves but not the 
other is due to that event having been corrupted by a second, 
or several more fold caustics moving into the source path as 
the caustic structure evolves. The drop in the proportion of 
total HMEs which are SHMEs in both of the time separated 
light curves relative to the proportion obtained using a static 
field is, as it must be, offset by an increase in the frequency 
of MHMEs. The proportion of cusp HMEs in the model light 
curve has been unchanged by the inclusion of proper motion. 
This is expected because these events are formed when the 
source moves through a high magnification region outside 
the parts of the source plane which are densely populated 
with caustics. 



4 CONCLUSION 

A new method has been developed to calculate the mi- 
crolensed light curve of an extended source. In realistic sim- 
ulations this method achieves a given accuracy for a light 
curve in less than one tenth of a percent of the time re- 
quired by the ray-shooting method. The gain in efficiency 
has been used to obtain extra resolution in realistic model 
light curves. This has allowed an investigation of different 
types of high magnification events to be made. 

Long term monitoring provides a method for measuring 
the physical parameters of the source and lens. As an exam- 
ple, the ratio of cusp to total HMEs is approximately twice 
as large for a source moving perpendicular to the shear than 
for one moving parallel to the shear. This result is indepen- 
dent of source size and could be used to obtain information 
on the direction of the source trajectory, a variable which 
is required in the determination of the transverse velocity. 
In addition, we note that SHMEs will be far more com- 
mon with respect to MHMEs for smaller source sizes. This 
result is approximately independent of trajectory direction 
and so could potentially be used to determine the value of 
the source size. 

Between 50% and 75% of the HMEs were found to be 



SHMEs in the static field model light curves. The introduc- 
tion of stellar proper motions into the microlensing model 
had the effect of lowering this fraction. However the drop was 
not significant - less than 20%. The relatively high fraction 
of SHMEs is a very exciting result since these events will 
allow the geometry of the continuum region of the source 
quasar to be measured from a single well sampled HME. 

Q2237+0305 has three images other than the one mod- 
elled here. Image B has parameters similar to those of image 
A, while images C and D are located in regions which have 
a higher optical depth. One direction for future work will 
be to repeat this analysis for model star fields which have 
optical depths corresponding to these other images. 



REFERENCES 

Chang, K., Refsdal, S., 1979, Nature, 282, 561 
Chang, K., 1984, Astron. Astroph., 132, 168 
Corrigan et al., 1991, Astron. J., 102, 34 
Fluke, C. J., Webster, R. L., 1999, M.N.R.A.S., 302, 68 
Foltz, C. B., Hewitt, P. C, Webster, R. L., Lewis, G. F., 1992, 
Ap. J., 386, L43 

Grieger, B., Kayser, R., and Refsdal, S., 1987, Astron. Astrophys., 
194, 54 

Haugan, S. V. H., 1996, IAU Sypmosium 173, p275. 

Irwin, M. J., Webster, R. L., Hewitt, P. C, Corrigan, R. T., Je- 

drzejewski, R. I., 1989, Astron. J., 98, 1989 
Katz, N., Balbus, S., Paczynski, B., 1986, Ap. J., 306, 2 
Kayser, R., Refsdal, S., Stabell, R., 1986, Astron. Astrophys., 166, 

36 

Kundic, T., Wambsganss, J., 1993, Ap. J., 404, 455 
Kundic, T., Witt, H. J., Chang, K., 1993, Ap. J., 409, 537 
Lewis, G. F., Miralda-Escude, J., Richardson, D. C, Wambs- 
ganss, J., 1993, MNRAS, 261, 647 
Lewis, G. F., Irwin, M. J., 1995, MNRAS, 276, 103 
Lewis, G. F., Irwin, M. J., 1996, MNRAS, 283, 225 
Paczynski, B., 1986, Ap. J, 301, 503 

Schmidt, W. R., Webster, R. L., Lewis, G. F. 1998, MNRAS, 295, 
488 

Schmidt, W. R., Wambsganss, J., 1998, Astron. Astrophys., sub- 
mitted 

Schneider, P., Ehlers, J., Falco, E.E., 1993, Gravitational lenses, 
Springer- Verlag, Germany 

Schramm, T., Kayser, R., Chang, K, Nieser, L., Refsdal, S., 1993, 
Astron. Astrophys., 268, 350 

Wambsganss, Kundic, T., 1995, Ap. J., 450, 19 

Wambsganss, J., Witt H., J., Schneider, P., 1992, Astron. Astro- 
phys., 258, 591 

Wambsganss, J., Paczynski, B., Katz K., 1990, Ap. J., 352, 407 
Wambsganss, J., Paczynski, B., Schneider, P., 1990, Ap. J., 358, 
L33 



© 1998 RAS, MNRAS 000, [I|-|h] 



10 J. S. B. Wyithe & R. L. Webster 

Wambsganss, J., Paczynski, B., 1991, Astron. J., 102, 864 
Wambsganss, 1992, Ap. J., 392, 424 
Witt, H. J., 1993, Ap. J., 403, 530 

Witt, H. J., Kayser R., Refsdal, S., 1993, Astron. Astrophys., 268, 
501 

Witt, H. J., Mao, S., 1994, Ap. J., 429, 66 



© 1998 RAS, MNRAS 000, 



